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We analyse a class of ID lattice models, known as Mfe models, which are characterised by an order- 
k clustering of spin-less fermions and by A/” = 2 lattice supersymmetry. Our main result is the 
identification of a class of (bulk or edge) defects, that are in one-to-one correspondence with so- 
called spin fields in a corresponding Zfc parafermion OFT. In the gapped regime, injecting such 
defects leads to ground state degeneracies that are protected by the supersymmetry. The defects, 
which are closely analogous to quasi-holes over the fermionic Read-Rezayi quantum Hall states, 
display characteristic fusion rules, which are of Ising type for fc = 2 and of Fibonacci type for fc = 3. 

PACS numbers: 05.30.-d, 71.10.Fd 


Introduction. - In the field of topological quantum 
computation (TQC) [1], a number of important lessons 
have been learned. The first is that non-Abelian statistics 
tend to be associated to a form of pairing or clustering 
in a quantum condensate. In 2D, p-wave pairing of spin¬ 
less fermions in a p + ip superconductor or Moore-Read 
(MR) quantum Hall state gives rise to non-Abelian statis¬ 
tics of Ising type, through the mechanism of Majorana 
bound states at the cores of half flux quantum vortices 
or quasi-holes [2, 3]. In the quantum Hall context, going 
beyond Ising anyons requires going beyond pairing, as in 
the Read-Rezayi (RRfe) [4] and NASSfc states [5]. The 
simplest examples beyond the MR state, the RR3 and 
NASS2 states, both give rise to Fibonacci anyons, which 
are universal for TQC (see, for example, [6]). More gen¬ 
erally, the anyons carried by the RRfc states are universal 
for fc = 3 and k > 5 [7, 8]. 

A second lesson learned is that a TQC-through- 
braiding protocol can be defined not just in 2D but 
also in a ID setting [9]. One starts from a T-shaped 
wire junction with non-Abelian defects at the wire ends 
and then runs a protocol of braiding, either in position 
space (by moving defects along the wires) or in parame¬ 
ter space. Again the prototypical example are Majorana 
bound states at the defect points. The underlying pair¬ 
ing is typically assumed to be extrinsic, meaning that it 
is induced (as in the Kitaev chain) through the proximity 
of a nearby superconductor. 

We here consider the question if one can construct ID 
lattice models with built-in, intrinsic, pairing or clus¬ 
tering properties and with defects binding Majorana 
or parafermion zero modes. We find that this goal is 
achieved by the lattice models M^ introduced in [10] and 
further analysed here. The definition of the M^ models 
involves a hard-wired /c-clustering constraint as well as 
Af = 2 supersymmetry. The order-A: clustering leads to 
Zfc parafermion degrees of freedom in the CFT describ¬ 
ing the Mfc models at criticality. In fact, the supercharge 
operator, which injects an extra particle into the sys¬ 
tem, contains the parafermion held ipi. A similar struc¬ 


ture is maintained if we move into a gapped phase. Our 
main result is the identihcation of a class of (bulk and 
edge) Mfc model defects that precisely correspond to the 
so-called spin helds Ci in the parafermion CFT. These 
defects are in many ways analogous to the quasi-holes 
over the RRfc quantum Hall states: they have fractional 
particle number and display characteristic non-Abelian 
fusion rules. The underlying mechanism is that of su¬ 
persymmetry protected order that is in essence of charge 
density wave (CDW) type. To turn this into supersym¬ 
metry protected topological order in the ID sense will 
require a non-local reformulation via a Jordan-Wigner 
type transformation. 

Mfc models. - The Mfc lattice models [10] describe 
spin-less fermions on a ID lattice, subject to the ‘order- 
k clustering’ constraint that at the most k particles can 
occupy consecutive sites. A supercharge Q+ is dehned as 

L 

f = l a,b 

where ^ is a fermionic creation operator which cre¬ 
ates a particle at lattice site j in such a way that a string 
of a particles is formed, with the newly created particle 
at position b. Choosing the such that = 0, 

we dehne a Af = 2 supersymmetric hamiltonian through 

H = {Q+,Q-}, (2) 

with Q~ = {Q'^y ■ This hamiltonian combines hopping 
terms with local potential and interaction terms. By con¬ 
struction, [id, Q+j = [H,Q~] = 0. All states in the spec¬ 
trum are doublets with [/, / + 1] particles, with the ex¬ 
ception of the supersymmetric groundstates at if = 0, 
which are annihilated by both Q^' and Q~. 

Possible choices for the have been studied in 

[10, 11]. Here we choose for the M2 model 

A[i,i],j = V2Xj, A[2,i]j = A[2,2].i = Aj, (3) 
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with the Xj staggered as ... lAlAl_ The factor y/2 

guarantees that the model is integrable [11] and, if A = 1, 
critical [ 10 ]. 

For the general Mfc model we choose parameters de¬ 
scribing a critical point perturbed by a specific, inte¬ 
grable, staggering [12]. The staggering, with lattice pe¬ 
riodicity k + 2, connects the critical regime with one of 
‘extreme staggering’ A <C 1, where the A[o^b]^j follow a 
simple pattern. For fc = 3, to lowest order in A, 
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with the dots indicating repetition modulo 5. We denote 
this as ... * *A * *..., with the ‘A’ indicating the central 
position in the staggering pattern. 

The Witten index for the M/j model with periodic 
boundary conditions (PBC) and with L = l{k + 2) sites 
is Wk = k + 1; indeed, for A > 0 the models have pre¬ 
cisely this number of supersymmetric groundstates, all 
ai E = 0 and filling v = k/{k + 2) [10, 13]. They are pro¬ 
tected against perturbations that commute with super- 
symmetry and do not affect the /c-clustering constraints. 
For open BC there are either zero or a single supersym¬ 
metric groundstate with E = 0, the latter for L = 0, — 1 
mod (fc -I- 2). We find, however, that in the presence of 
suitable boundary or bulk defects, the open systems have 
states with energies that are exponentially suppressed, 
E oc with Li characteristic distances among de¬ 

fects and boundaries. The exponential degeneracies are 
protected by supersymmetry. 

M2 model. - We now zoom in on the M2 model on 
an open chain. At criticality (A = 1), the finite size 
spectra can be matched with those of the 2 nd minimal 
model of A/” = 2 superconformal field theory (CFT), of 
central charge c = |. The match can be made with 
the help of numerical spectra (we analysed open chain 
spectra up to length L = 25, fig. 1) and are similar to 
the results of [14] for the Mi model. The relevant CFT 
modules are Vm, ipVm with m G Z + ^ and aVm with 
m G Z. Here the Vm are charge m vertex operators for 
a c = 1 scalar field and the ip, a arise from the c = | 
Ising CFT factor. States in Vm have an even number 
of ^-modes while those in ipVm contain an odd number. 
The supercharges are the (Ramond sector) zero modes of 
the supercurrents ipV± 2 (z). 

Up to an overall 1/L scaling, the lattice model energies 
correspond to Ecft = Lq — The lowest energies are 
(4m^ — 1)/16 for Vm, (4m^ -I- 7)/16 for ipVm and m^/4 
for aVm- The lowest-energy states in U_j_i and crVb are 
supersymmetry singlets with Eqft = 0, all other states 


have Ecft > 0 and pair up into doublets. The critical 
M 2 spectrum with open BC is easily described: with 
m = 2/ — L — i, one finds the CFT modules Vm for 
/ even and ipVm for / odd (fig. 1 ). 

Reducing A below 1 sends the theory off criticality, 
with RG flow leading to the supersymmetric sine-Gordon 
(ssG) theory at coupling = Stt [11]. The M 2 model off- 
critical finite size spectra can be analysed in terms of ssG 
bulk S'-matrices and boundary reflection matrices [15]. 
The ssG theory holds important clues for the topological 
aspects of the M 2 model degeneracies [15-17]. 

Rather than following the RG flow, we will here con¬ 
sider the limit A ^ 1 (‘extreme staggering’), where the 
M 2 eigenstates approach a simple factorized form. This is 
analogous to a special tuning in the Kitaev chain, which 
leads to perfectly decoupled Majorana edge states [9]. 
This simple setting enables us to demonstrate how dif¬ 
ferent BG result in exponential ground state degeneracy 
beyond this idealised limit. The A ^ 1 limit is also sim¬ 
ilar to the thin-torus limit of the MR and RRfe quantum 
Hall states [18-20]. Indeed, the systematics of the fusion 
channel degeneracies is highly analogous between the two 
settings. 

For A = 0 and PBG, the M 2 groundstates are 

]-) = ...0(T.)0(T..., ]+) = ...l.)0(-l-)0..., 

10) = ...1010101..., (5) 

where (T-) = 110 -I- 011 and ]—) and J-l-) are related by 
a shift over two lattice sites. For open BG, there is at 
the most a single E = 0 groundstate for given particle 
number /. For L = 4^ — 1, staggering lA... Al, / = 21, 

l+)o,o = [(-l.)0(-l-)---(-l-)] , (6) 

where ‘ 0 , 0 ’ refers to open/open BG. For A > 0 this state 
remains at E = 0, where it is protected by the Witten 
index, W = 1, and it is separated from all other states 
by a gap that remains finite as long as A < 1. 

Boundary defects. - To steer into a case with expo¬ 
nentially degenerate groundstates at given particle num¬ 
ber /, we need to enforce a defect at both boundaries 
that allows all three A = 0 PBG groundstates to connect 
to the edge at zero energy cost. For this we impose the 
constraint that the two sites adjacent to a boundary can¬ 
not both be occupied by a particle [21]. With this BC 
(which we call of ‘cr-type’ and denote by a bracket .. .]o-), 
all three A = 0 PBC vacua can connect to the boundary 
at zero energy cost. This gives, for L = 41 + 1, Al... lA, 

I-).,. = .[O(T-)... 0 ],,, = ,.[ 100 (- 1 -)... 001 ],,, 

|0)..„ = .[1010...0101]„. (7) 

For A > 0 the state \—)a,a at f = 21 remains at U = 0 
while \-\-)a,a at f = 21 and [Ojo-.cr at / = 21 -|- 1 pair into 
a doublet of energy 6E{X) > 0. The energy difference 










3 


0 


X V'^-| I }5 

}3 


X 


^ }4 


X v_3 X IpVl 


Vs 


0 


aV.3 ^ }i 


n}4 

X } 2 X } 2 

X (jVIi X aVi 


X o-Va 


0 


X 
X 

X 

X w-i i 

X V_s ^ 
2 

X 


W- 

V_i 


X 

X IpVs 

,1 X V^3 ' 


X ^Vt 
X Vr^ 


-7/2 


-5/2 


7/2 


5/2 "3 -3-1 1 3 371 

(b) BC a jo-pen 

FIG. 1. Numerical M 2 spectra with L = 25, / = 11,12,13,14 up to i? = 4. The labels specify the corresponding CFT modules. 


-3/2 1/2 

(a) BC openjopen 


-1/2 3/2 

(c) BC aja 


between originates from the boundary, which im¬ 

plies that it involves a power of A that scales with the 
length of the system. We checked numerically that 

8E[\) oc (8) 

which gives exponential splitting 6E cx e~°‘^ as long as 
A < 1. We checked that this behaviour is robust against 
perturbations obtained by deforming some of the param¬ 
eters \a,b],j in eq- (l)j provided we do not break the 
Af = 2 supersymmetry. We thus identify supersymme¬ 
try as the agent protecting the exponential degeneracy 
of our ‘qubit’ |±)ct.(t- We remark that many ‘natural’ 
perturbations of the M 2 model do break supersymmetry. 
An example is the local fermion density operator pi at 
site i = Ap + 1. AtA = 0 the expected value is (pi) = 0 
in the state |—)c 7 ,cr while (pi) = 1 in the state |+)c 7 ,cr- 

Following the ‘qubit’ states all the way to the CFT 
point, A = 1, we find (fig. 2) 

\-)a,a l+)o',cr (9) 

where I-) denotes the lowest weight state of the corre¬ 
sponding CFT module. At the CFT point the boundary 
Majorana modes have delocalised and the energy split¬ 
ting is of order 1 /L. 

The degeneracy of the ‘qubit’ |±)ct.(t can be traced to 
the fusion channel degeneracy aa = 1 ip oi the Ising 
spin-field a in the underlying CFT. Through the qH-CFT 
connection [3], this choice of fusion channel carries over to 
the fusion product of two quasi-holes over the MR state. 
We consider the MR state in spherical geometry, which 
we view as an open ‘tube’ capped by specific boundary 
conditions at the two poles. For N = 21 particles, the 
MR groundstate has the following thin-torus form 

MR,iV = 2/: [11001100... 110011] . (10) 

The analogous groundstate of the M 2 model is precisely 
the state |+)o.o in eq. ( 6 ). The simplest case with two¬ 
fold fusion channel degeneracy is that of the MR states 



FIG. 2. Flow between ‘extreme staggering’ (left, A = 0) and 
critical (right, A = 1) limits of the M 2 model with i = 17 
sites, / = 8 particles, staggering A1. .. lA and aja BC. The 
lower two states constitute the ‘qubit’ eq. (9). 


with 71 = 4 quasi-holes. The general counting formula for 
71 quasi-holes reads [ 22 ] 


E 

F=N mod 2 


N-F 

2 


+ n 


n 



( 11 ) 


Here the first binomial counts orbital degeneracies of the 
77 quasi-holes, while the second, together with the sum 
over F, pertains to the fusion channel degeneracy. We 
fix the orbital degeneracy by selecting the states with two 
quasi-holes at both the north and the south poles. 


F = 0: [01100... 110], F = 2: [10011... 001]. (12) 

These ‘MR qubit’ states correspond to the ‘M 2 qubit’ 
states |±)ct.(t of eq. (7). To understand this we have to 
compare the open M 2 chains with the MR states not on 
the sphere but on the cylinder. On a cylinder with vacua 
‘1100’ at far left and right, we can extend the F = 0 state 
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as 

... 1100^^[01100... 110]<,^0011..., (13) 

where a a denote the two quasi-holes at the boundaries. 
We can move one of the quasi-holes out from each of the 
boundaries to get 

... 1100^1010... 10^[01100... 110]<,0101... 01^0011.... 

. 

This corresponds to the situation that we have in the 
M 2 model, where cr-type BC arise from the presence of 
a single a quantum at a boundary. This interpretation 
is confirmed by the CFT content of the open chain finite 
size spectra at A = 1 (fig. 1). For open/cr BC, we find 
all sectors aVm with m shifted to m = 2f — L. For a/a 
BC, with non-Abelian defects at both ends, we find, at 
m = 2f — L + ^, both the Vm and i/Vm modules, in 
accordance with the fusion rule aa = 1 + 1 /. The states 
in eq. (9) are a particular example. 

Bulk defects and quantum register. - At extreme 
staggering, kinks connecting any two of the A = 0 vacua 
come at finite energy cost. For example, kinks/anti-kinks 
connecting |0) and |-|-), written as 

K K 

...lOlOl'o 0 (-l-)..., ...lOlOlTo(-l-)..., (15) 

have E = 1] the same holds true for all kink types (a, b) = 
(0, ±),(±,0). Note that kinks/anti-kinks are connected 
by A/” = 2 supersymmetry, 

g+ : Ka,b ^ iFa.fe. ( 16 ) 

Multi-kink/anti-kink states can be counted through for¬ 
mulas similar to those for the MR state, see eq. (11), 
for all choices of open chain BC. Followed through to the 
CFT limit, these counting formulas provide novel expres¬ 
sions for characters of the CFT [15]. 

We now define cr-type bulk defects. These will allow 
some of the bulk kink states to exist at zero energy (for 
A = 0). In the example of eq. (15) this can be done by 
excluding the configuration ‘11’ at the kink location: this 
eliminates the anti-kink and turns the kink into an if = 0 
state! To treat the ± vacua on equal footing, we repeat 
the same constraint two steps to the right, and define 

cr:...Al Al..., cr':...lA lA.... (17) 

I_I I_I I_I I_I ' ’ 

no 11 no 11 no 11 no 11 

This definition holds for general A and represents a con¬ 
straint in the Hilbert space of states. At A = 0, a defect 
a can connect a vacuum |0) coming in from the left to 
|-|-), |0) or |—) extending to the right, and opposite for 
cr', all at zero energy. 

We can now conceive a ‘quantum register’ by taking an 
open chain, length L = EL + 1, staggering type Al... lA, 


CT-type BC at both ends, and injecting a sequence of 2n 
well-separated defects in the order a'a ... a'a. This leads 
to 3”+^ degenerate groundstates at A = 0, with 2”+^ of 
them having the minimal particle number f = 21 — n. 
These 2”+^ states form an Ising anyon ‘quantum regis¬ 
ter’. At A > 0, a unique E = 0 groundstate at f = 21 —n 
remains, while all other states pair up into doublets at 
energies of order 

M 3 model. - Turning to /c = 3 we identify the follow¬ 
ing four PBC groundstates at extreme staggering 

|1) = ...11100..., |i) = ...(-l---)..., 

|0) = ...01011..., |-i) = ...0(Tl-)..., (18) 

with ^ indicating the position of ‘A’ in the staggering 
pattern eq. (4), (-1 • • •) = 01101 - OHIO -f 11001 - 11010 
and (T1-) = 1110 — 0111. The energies of the kinks 
ATa.b are nia^b = 2|a — b\. This is in agreement with 
the kink masses in the Af = 2 supersymmetric massive 
integrable QFT with Chebyshev superpotential W{X) = 
— P'^X^ + j3X [23-26], which we expect to result 
from the RG flow set by the staggering perturbation. 

On open chains, we define cti type BC by excluding 
‘ 111 ’ near a given edge, and a 2 type by excluding ‘ 11 ’. 
At the level of the open-chain CFT spectra, a ai type BC 
precisely corresponds to the Z 3 parafermion spin field 
ai'. upon changing BC, the various CFT sectors shift 
according to the fusion products with these two spin 
fields. The k = 3 CFT contains, besides a free boson, 
Z 3 parafermions ipi ^2 and parafermion spin fields ai, e. 
The supercharge Q"*" is the zero-mode of the super cur¬ 
rent il;iV^{z). For open/open BC, the M 3 model realises 
the sectors, with m= §/ — fA — |, 

{VW,V’lKn,^ 2 K„} (19) 

for fc = 0,1, 2 with k = f mod 3. For open/cri BC this 
becomes, with m = 

{aiVm,eVm,cr2Vm}, ( 20 ) 

in accordance with the parafermion fusion rules aitpi = s 
and criV ’2 = o' 2 - For open/CT 2 BC, with m = — ^L+ A, 

{cr2Vm,criVm,£Vm}, ( 21 ) 

in agreement with a2'4’i = and a 2 'i /’2 = £■ The super- 
symmetric groundstates are |cri_ 2 l±i) and iH-ta). 

Putting ai type BC on both ends, the open-chain CFT 
sectors follow the fusion rules ctiCTi = i/i + a 2 , aia 2 = 
1 + e and a 2 a 2 = 1/2 + <^ 1 - As for fc = 2, these fusion 
rules lead to exponential groundstate degeneracies in the 
extreme staggering limit. A characteristic case wonld be 
L = 15 sites, (T 2 -type BC on both ends, and staggering 
positioned as **A .... Here the lowest CFT states are two 
supersymmetry doublets (at / = 8,9) at CFT energies 
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FIG. 3. Fusion rules of bulk defects in the M 3 model at ex¬ 
treme staggering. 


^'CFT = l/5j ^'CFT =4/5, while the A = 0 M 3 model 
has four E = Q vacua. The states connect according to 

\\)<T 2 ,a 2 ^ IV'i 4T5/4), |0)cr2,(72 ^ IV’ 24 + 5 / 4 )) 
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We define bulk defects as 
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Fig. 3 depicts the corresponding zero-energy fusion rules. 
They determine the size and structure of the quantum 
register opened up by an alternating sequence of defects 
(Ti, a'y In all cases the states at maximal particle number 
are made up entirely of vacua |0) and |1). Restricting 
the fusion rules to |0), |1) (drawn lines) gives Fibonacci 
number degeneracies. The Af = 2 supersymmetry acts 
within the register, with Q~ mapping the |0), |1) into 
linear combinations of |4), | — 4). In the example of L = 
51 sites, cF 2 l<J 2 BC, staggering type ** A ... and n cr 2 /cr 2 - 
type defects, the degeneracy at / = 31 is Fibo„+ 3 , with 
Fiboj = 1,1,2,3,— Eq. (22) is the case n = 0. For 
n = 1 one finds 3,4,1 states at f = 31,... ,31 — 2, n = 2 
gives 5, 9, 5, 1 states at f = 31,... ,31 — 3, etc. 

For the general M^, model, defects eliminating fc-|-l —j 
consecutive T’s will correspond to the parafermion 
spin fields cr^, j = 1,..., fc — 1 [15]. 

Conclusions. - We have demonstrated that (bound¬ 
ary and bulk) defects in the models off-criticality lead 
to quantum registers that are protected by supersymme¬ 
try. It will be interesting to explore ways to manipulate 
these registers, so as to act on the quantum informa¬ 
tion that can be stored in the supersymmetry protected 
ground states. One idea is that of a ID braid protocol (as 
in [9]); one expects that this will result in the braid ma¬ 
trices as they are known for the corresponding 2D (RRfc) 
topological phases. It will also be interesting to see if a 
dual formulation of the M^, models, with ID topological 
order rather than order of CDW type, can be obtained. 
One expects that operators that preserve supersymmetry 


in the models correspond to operators that are local 
in the dual variables of the topological phase, and are 
thus unable to split the exponential degeneracies. 

Many extensions of the ideas presented here are feasi¬ 
ble. In the M 2 model, alternative bulk defects, based on 
a ‘no 0’ rather than a ‘no 11’ condition, lead to Fibonacci 
number degeneracies. The Mi model on a square ladder 
is in many ways similar to the M 2 model [27] and it is 
natural to introduce cr-type defects in the 2D Mi model 
on the octagon-square and square lattices [28, 29]. In 
all these cases, we expect non-trivial fusion relations to 
emerge. 
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